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Abstract 

We present experimental results for pattern formation in Rayleigh-Benard 



> 

o 
o 

(•^ . convection of a fluid with a Prandtl number a ~ 4. We find that the spiral- 

ly . 
^N . defect-chaos (SDC) attractor which exists for cr ~ 1 has become unstable 



Gradually increasing the temperature difference AT from below to well above 
its critical value ATc no longer leads to SDC. A sudden jump of AT from below 
to above ATc causes convection to grow from thermal fluctuations and does 



^ . yield SDC. However, the SDC is a transient; it coarsens and forms a single 

cell-filling spiral which then drifts toward the cell wall and disappears. 
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Rayleigh-Benard convection (RBC) occurs in a shallow horizontal layer of a fluid heated 
from below when the temperature difference AT across the layer exceeds a critical value AT^. 
The velocity fields which form above ATc have long been used as an experimental testing 
ground for theories of pattern formation in nonlinear systems. [^ The phenomena which 
occur depend upon the Prandtl number a = zz/k, the dimensionless ratio of the kinematic 
viscosity u to the thermal diffusivity k. A fascinating recent discovery in RBC is that of 
spiral-defect chaos (SDC) for cr ~ 1. |^ It consists of patches of rotating convection-roll 
spirals and of other defects, and occurs for e = AT/ ATc — 1 > e^ > 0. The spirals are 
coherent structures with a typical diameter of a few roll wavelengths. They can be right- 
or left-handed, single- or multi-armed, and appear and disappear irregularly in time and 
space. SDC has since been found in numerical solutions of model equations [^ and of the 
Navier-Stokes equations (NSE) Q. The band of wavenumbers for SDC 0,^ lies well within 
the range for which straight rolls are also stable p[ in the laterally infinite system. Thus the 
solutions of the equations of motion of RBC exhibit bistability of time-independent straight 
rolls and of SDC, with each of the two having distinct attractor basins. For cr ^ 1 the generic 
initial conditions of real RBC experiments lie in the attractor basin of SDC. 0,^ Here we 
report that this situation is changed dramatically when a is increased. For a = 4.0 ± 0.2, 
we do not find any spirals when increasing e in small steps from below zero to large positive 
values. Thus, this experimental path does not cross an attractor basin of SDC. However, 
when e is increased suddenly from below zero to values of 0{1), SDC evolves in the sample 
interior from random noise 0. Thus, random initial conditions appear to remain within the 
attractor basin of SDC even for o" ~ 4. But this is illusory since the SDC turns out to be 
a transient rather than a stable attractor. At very long times, the spirals coarsen, forming 
larger and larger spirals until only one giant spiral fills nearly the entire cell. At even longer 
times the giant spiral is not stable either; its head slowly drifts toward the sidewall and 
disappears from the cell, leaving a state similar to those obtained by increasing e gradually. 
Our observations differ from those of Assenheimer and Steinberg, [^^ who reported SDC 
in a sample of SFg near its critical point with Prandtl numbers similar to and even larger 



than ours. 

A detailed description of the apparatus was given in Ref. [|T0|. The convection cell had 
a 0.95 cm thick sapphire top plate and a diamond-machined aluminum bottom plate. A 
film heater glued to the lower side of the bottom plate provided the heat current. The top 
of the sapphire plate was in contact with a temperature-regulated circulating water bath. 
The temperatures for both bottom and top plates were regulated to better than 1 mK. The 
circular sidewall was made of high-density polyethylene. It had an inner diameter of 8.89 
cm, and was sealed to the top and bottom plates with ethylene-propylene 0-rings. High- 
purity liquid-chromatography grade acetone was dried by using SA molecular sieves followed 
by distillation, and then injected into the cell with a syringe without significant exposure 
to the atmosphere. Two cells were used. The cell height was d = 0.80 ± 0.01 (1.0 ± 0.01) 
mm, corresponding to an aspect ratio (radius/height) F ^^ 55 (45). The calculated value 
of ATc was 8.56 (4.50) °C. At onset, the coefficient Q ||ll[ describing the departure from 
the Boussinesq approximation was 0.37 (0.20). The vertical thermal diffusion time ty was 
6.8 (10.5) s. The onset of convection was determined from measurements of the Nusselt 
number Af (the ratio of the effective thermal conductivity to that in the conduction state). 
Quasi-static results for A^ over the entire range of the experiment are shown in Fig. ^ Images 
of convection patterns were taken using shadowgraph apparatus described previously ||12 . 



In one set of experiments, e was varied slowly from below zero to above 5 by changing it in 
steps of about 0.1. The system was allowed to equilibrate for at least one horizontal diffusion 
time V^ty (about 6 hours) after each step. This procedure corresponds to an effective ramp 
rate f] = de/dt < 5.4 x 10^^ (here t is measured in units of t^). It yielded images like 
those shown in Fig. |^. For e ^ 1, the patterns consisted essentially of straight rolls, except 
for typically one or two grain boundaries and perhaps one focus singularity adjacent to the 
sidewall. A representative example is shown in Fig. ^, which is for e = 0.5. The wide 
range of e over which relatively straight rolls occur is in contrast to the a ~ 1 case, where 
significant roll curvature occurs already for e > 0.1 |]T3|,|T1|. Figures ^, c, and d are for e = 



2.0, 3.7, and 4.5. They show that the rolls do acquire significant curvature as e is increased. 



This curvature is associated with a tendency for the roll axes to terminate perpendicular 
to the wall. It leads to the formation of additional wall foci and to a reduction of the roll 
wavelength in the cell center. For a ~ 1, this phenomenon occurred already at e values 
which were about an order of magnitude smaller, and yielded a time- dependent state with 
repeated formation of defects in the cell interior via the skewed- varicose mechanism. |TB|JT^ 



In the present case, the images were essentially independent of time at constant e. No defect 
formation occurred because even the significantly enhanced wavenumber range encountered 
in the cell interior did not result in any crossing of the skewed- varicose instability-boundary. 
We illustrate this latter point in Fig. |], which at small wavenumbers gives the Eckhaus 
and at large the skewed-varicose instability as a solid line. The solid circles are the mean 
wavenumbers k from the present run. They were obtained from the structure functions S{k) 
(squares of the moduli of the Fourier transforms) of the images as described elsewhere 0. 
The horizontal bars represent the widths of the structure functions, obtained from the second 
moments of S{k) about k. It is apparent from these data that even the large- and small-/c 
tails of S{k) are primarily well within the range of stable straight rolls of the infinite system. 

When e was increased further, the formation of cellular flow began near one of the wall 
foci, as shown in Fig. ^ for e = 5.0. During many horizontal diffusion times these flow 
cells multiplied and spread over a wider part of the sample and evolved into a pattern of 
coexisting regions of hexagons with either upfiow or downflow at their centers. These regions 
occurred together with curved rolls and wall foci in other parts of the sample. This is shown 
in Fig. ^, which is also for e = 5.0 but GOOOt^ after e was last increased. This pattern is 
similar to ones observed by Assenheimer and Steinberg (AS) with SFg close to its gas- 
liquid critical point in a cell with F = 80 at e = 3.4 and a = 4.5. However, the threshold 
for and mechanism of hexagon generation seems different in our case, since AS report that 
their hexagons appear already near e = 2.8, and via a core instability of spirals or targets 
which do not exist in our system. Apparently the threshold for hexagon generation depends 
on the nature of the previously existing structure. 

Of greatest relevance to the present issue is that none of the images obtained in runs like 



the one described above revealed any spirals. This is so also for the larger F = 55. In contrast 
to this, for cr ~ 1 and similar values of F [0,^] the onset of SDC was found at e^ ~ 0.55, 
and SDC persisted to above e ~ 3 where the oscillatory instability was encountered at that 
a value. The total absence of spirals in our experiment also differs from the results reported 
by AS, who for a = 4.3 found spirals and targets (zero-armed spirals) for e ^ 0.84. 

In a second set of experiments we increased e suddenly from eo = — 1 to ei > 0. Care 
was taken to increase e as quickly as possible to the final value in a smooth, monotonic 
manner without overshooting. This procedure yielded small-amplitude random convective 
flows in the cell center after only a few ty. Examples are shown in Fig. ^ja and e. The 
dynamic sidewall forcing during the transient typically led to a small number of concentric 
rolls surrounding the random flow field in the center. Stabilizing at various values of ei ^ 1, 
we observed a rapid growth (typically within SOt^) of targets and spirals from the disordered 
flow field, as shown in Fig. ^d and f. Simultaneously the mean wavevector decreased and 
after a relatively short time settled down at the points shown by open circles in Fig. ^. Once 
the pattern of spirals, targets, and other defects was generated, the complicated dynamics of 
SDC known from the experiments at smaller a came into play. For smaller ei (say ei ~ 0.5), 
no spirals or targets developed and a pattern like that of Fig. ^ evolved after a relatively 
short time. 

In contrast to the experiments with cr ~ 1, the SDC state which for ei ^ 1 grew from the 
random initial conditions was not stable. Instead of persistent SDC dynamics, we observed 
a coarsening of the patterns over time. This is illustrated by Figs. ^ and g. The spirals 
grew in size until after a few horizontal diffusion times there was left primarily one giant, 
nearly cell-filling, slowly rotating spiral as shown in Figs. ^ and h. 

Although the giant spirals lived for several horizontal diffusion times, experiments over 
even longer times showed that they were not really stable. This is illustrated in Fig. |^. At 
7500t„ after a jump in e to ei = 1.6, one giant spiral is positioned with its core centered in 
the cell. As shown by the other images, this core location is not stable; the core gradually 
drifts towards the wall while the radius of the spiral decreases so that it roughly equals 



the distance of the core from the walL During this process, no new spirals or targets are 
born. At a time shortly after that of image d, the core left the cell altogether and a state 
of straight/curved rolls and wall foci similar to the ones obtained by the slow increase of e 
(Fig. 1^) remained. 

For cr ~ 4 and F = 45 or 55, we found that SDC can be attained from the random flow 
fields generated by a jump in e, but not by quasi-static changes of e. This is in contrast 
to cr < 1 where SDC is difficult to avoid for F > 30 p|,P no matter how e is varied. We 
also demonstrated by experiments extending over very long times that SDC, once formed, 
is unstable in our samples with a ~ 4. It coarsens and evolves into a single cell-filling 
spiral, whose core then drifts from the cell center towards the sidewall where it disappears. 
Our results differ from those of Assenheimer and Steinberg P,^, who found SDC without 
coarsening over the wide range 2.8 ^ cr < 28, with the spirals gradually being replaced by 
targets as a increases. The reason for this difference is unclear. However, we note several 
differences between the two experiments. In one run with cr = 4.5 described by AS in some 
detail [^] they used a cell for which AT^, was two orders of magnitude smaller than ours. 
Thus, assuming that their temperature control was similar to ours, their experimental noise 
(spatial as well as temporal) in e presumably was two orders of magnitude larger than ours. 
Furthermore extremely slow ramp rates /? ~ 5 x 10~^ like ours would be more difficult to 
attain when AT^ is very small. Also, the total time intervals over which they observed SDC 
at constant e were only (9(F^r^) [jl5[ and are at least an order of magnitude smaller than 
ours (~ 18F^t^). Finally we note that AS used a sample with F = 80 which is a factor of 1.5 
larger than one of ours. One might conjecture that larger F (= 80, the AS value) stablizes 
SDC for a ~ 4 whereas F = 55 (our larger value) destabalizes it. This would require a 
strong a dependence of €^(0", F) because SDC is readily attained for instance for F = 30 and 
cr = 1. If this trend in a dependence continues for o" > 4, AS should not have observed 
SDC in their cell with F = 80 at their much higher Prandtl numbers, up to 20 or larger. 
We conclude that the differences between the two experiments remain unresolved, and that 
additional measurements, particularly as a function of F, are required to shed light on this 



problem. Unfortunately it is difficult to reach larger F values using acetone because ATc 
becomes very large if d is reduced. 

Finally we remark that it would be interesting to search for the coarsening and instability 
of SDC in solutions of the Navier-Stokes equations. 

We are grateful to Werner Pesch for permitting us to use his program and for his advice 
during the calculation of the stability boundary in Fig. ^. This work was supported by the 
National Science Foundation through Grant DMR94-19168. 
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FIGURES 
FIG. 1. The Nusselt number Af vs. e for F = 45 obtained while increasing e quasistatically. 

For this run, the top temperature was 16.5°C and ATc = 4.50 °C. 

FIG. 2. Representative images obtained during the experimental run of Fig. |l| where e was 
increased gradually, a) e = 0.5; b) e = 2.0; c) e = 3.7; d) e = 4.5. Up to this value of e, "rolls" 
(with curvature, wall foci, etc.) are stable and time independent. Images e) and f) are for e = 5.0, 
lOOOti, and 6000i^ after increasing e. They show the co-existence of hexagon domains with upflow 
and downflow at the cell centers and with rolls and wall foci. No spirals or targets were encountered. 

FIG. 3. The solid circles are the mean wavevectors obtained in the run of Fig. || by increasing 
e quasistatically for F = 45. The open circles are wavevectors obtained from the SDC state before 
the coarsening process had progressed significantly. The horizontal bars through the data points 
correspond to the width of the structure function S{k). The plus and cross are the wavevectors 
reported in Ref. |^ for rolls and hexagons respectively. The solid line is the stability boundary for 
straight rolls in the laterally infinite system. 

FIG. 4. The time evolution of the patterns after changing e suddenly from eo = — 1 to ei > 0. 
Images a) to d) are for F = 45 and ei = 2.0 with a top temperature of 16.5 °C and ATc = 4.50 
°C. In units of t^ = 10.5 s, the elapsed times since the jump in e are a) t = 2, b) t = 400, c) t 
= 2000, and d) t = 8000. Images e to h are for F = 55 and ei = 1.0 with a top temperature of 
19.0 °C and ATc = 8.56 °C. In units of t^ = 6.8 s, the elapsed times since the jump in e are e) t 
= 2, f) t = 600, g) t = 6000, and h) t = 10000. In both examples the pattern in the cell interior 
developed from random convective flows, and the spiral/target structures coarsened and healed to 
a single spiral and some defects along the wall over a few horizontal diffusion times F^t„. 

FIG. 5. The long-time evolution of the patterns after changing e suddenly from eo = —1 to 
ei = 1.6 for F = 45. For this run, the top temperature was 18.0 °C, and ATc = 4.53 °C. The 
elapsed times (in units of ty) since the jump in e are a) t = 7500, b) t = 19000, c) t = 30000, and 
d) t = 38000. 
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